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Information of the phase structure of many flavor SU(3) gauge theory is of great interest for 
finding a theory which dynamicaUy breaks the electro- weak symmetry. We study the SU(3) gauge 
theory with fermions for Nf — 12 and 16 in fundamental representation. Both of them, through 
perturbation theory, reside in the conformal phase. We try to determine the phase of each theory 
non-perturbatively with lattice simulation and to find the characteristic behavior of the physical 
quantities in the phase. HISQ type staggered fermions are used to reduce the discretization error 
which could compromise the behavior of the physical quantity to determine the phase structure at 
non-zero lattice spacings. Spectral quantities such as bound state masses of meson channel and 
meson decay constants are investigated with careful finite volume analysis. Our data favor the 
conformal over chiral symmetry breaking scenario for both Nf = 12 and 16. 
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1. Introduction 

The technicolor model is one of attractive candidates of the physics beyond the standard model 
for the dynamical origin of the electroweak symmetry breaking. In particular the gauge theory with 
walking behavior near the edge of the conformal window has been considered as a realistic new 
physics model (walking technicolor). Therefore it is important to understand the phase structure 
and the low-energy dynamics of these theories for the phenomenology to be tested by the on-going 
LHC experiment. A concrete example of such models could be the SU(3) gauge theory with 12- 
flavor of massless fermions in the fundamental representation. The main purpose of this study is to 
understand the dynamics of this gauge theory by using the technique of the numerical simulation 
developed in lattice QCD, and ultimately to test whether this theory is a candidate of the walking 
technicolor (see [|l|] for more discussion and references). 

A number of lattice works studied this theory, such as the ones investigating the phase diagram 
of Wilson fermions [^, ^ , hadron spectrum [^, ^ ^ , and the running coupling constant ^, ^ ^] . 
However, whether this theory is in conformal or chiral broken phase is controversial at present. 

We investigate the 12-flavor SU(3) gauge theory using a variant of the highly improved stag- 



gered quark (HISQ) action [11], eventually with multiple lattice spacings, aiming to shed light 
on the controversy. We present preliminary results of pseudoscalar meson mass and decay con- 
stant, through which we attempt to determine the phase of this theory, as well as the characteristic 
quantity associated with the phase. 

Along this line we also simulate the 16-flavor SU(3) gauge theory with the same improved 
action. This theory resides deep in conformal window according to the perturbation theory. If this 
is the case, this theory could be used as a reference for the signal of conformality. Prehminary 
results for this study is presented here. 

Another series of simulations with this set up for Nf < 12 have been performed to test the 
walking behavior of Nf = 8, and to have a reference for the signal of the real QCD like hadron 
phase at Nf = 4. The results are reported in these proceedings 

2. Simulation setup and method 



We adopt the highly improved staggered quark (HISQ) action Jl IQ, but without the tadpole 
improvement and the mass correction term for heavy quarks. This action further suppresses the 
taste breaking in the QCD simulations compared to the Asqtad or stout-smeared staggered fermions 



[12]. We use this action for the many flavor simulation to minimize the discretization effects which 
potentially compromise the behavior of the physical quantities at non-zero lattice spacings. 

Gauge configurations are generated through HMC algorithm with various parameter sets for 
the fermion mass my, volume and the bare coupling j8 = for Nf = 12 and 16. We measure 
the mass and decay constant of the lowest state of the pseudoscalar channel (pion), and fj^, 
respectively, which are tested against the chiral symmetry breaking or conformal scenario. The 
tree-level chiral perturbation theory (ChPT) is used for the former. The conformal hyperscaling 



[ ]13| ] and and its finite-size version is used for the latter. 

The finite-size hyperscaling [|T^] is derived in the conformal theory deformed by the small 



fermion mass and also the large, finite volume L^. This scaling explains the fermion mass and 
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volume dependence of the physical quantities, such as the hadron mass mh, which is governed by 

the mass anomalous dimension at the infrared fixed point 7*. According to this scaling, a physical 

1 

quantity is described by a function /(a) of the scaling variable x = Lm^^' and L, 

mH=f{x)/L. (2.1) 
Since the finite-size hyperscaling is an extension of the hyperscaling to a finite volume, this scaling 



should reproduce the original hyperscaling [ |130 in the infinite volume limit, 



mH = CHmy'\ (2.2) 

Therefore when the L is large enough, the function should become 

/(x) = CO + C//X. (2.3) 

We use this fit form in the following analysis, and check if our data present the scaling behavior or 
not. 

3. Results of 12 

After a pilot study with several bare gauge couplings j3 = for the A^^ = 12, we determined 
to carry out the production run for j3 = 3.7. We generate the gauge ensembles at three volumes, 
X T, {L,T) = (12,24), (18,24), (24,32) \ and the various fermion masses, 0.04 < m/ < 0.2. 
We accumulate typically 500-1000 trajectories at each parameter set of L and nif, and measure the 
ntjc and fjc on these configurations. 

3.1 ChPT vs. hyperscaling fit 

Figure [1] shows the result of the niji at each L and m/. In this analysis we use the result obtained 
from only the largest volume, L = 24, because in this volume we expect that finite volume effects 
of the is negligible in our fermion mass range. 

We first attempt to analyze pion mass with a fit form motivated by the chiral perturbation 
theory (ChPT) given by 

ml = co + cimf + C2mj. (3.1) 

As the nij,; should vanish at the chiral limit, cq must be zero within the error if the chiral symmetry 
is broken, thus, ChPT describes our data. The values of the ;^^/d.o.f. are 4.1 and 28 for the fit form 
with and without the cq in eq. (^?T|), respectively. The former gives better ;^^/d.o.f, while the cq is 
non-zero cq = —0.0287(26). On the other hand, the hyperscaling fit eq. ( |2.2| ), gives ;^^/d.o.f. = 3.6 
and 7, ~ 0.46. This means that the hyperscaling, where the fit gives the smallest ;^^/d.o.f. with 
rrijc = at the chiral limit, is favored by our data. 

3.2 Analysis with finite-size hyperscaling 

1 

Figure ^ shows the Lmj^ at three lattice sizes as a function of the scaling variable x = Lm'j-^^ , 
with 7 being different in each panel. In 7 = 0, the data are scattered, but as the 7 increases, the 

'ideally, the ratio L/T should be fixed, which is not the case for L = 12. However, we find the spectrum is not 
sensitive to L/T when mj^L ^ 1, and we only use the data in such a region in this analysis. 
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Figure 1: in the 12 flavors as a function of the bare fermion mass nif. The different symbols denote 
the data at the different volumes. The solid and dashed lines represent the ChPT fit with and without the 



constant term cq in eq. (3.1 ), respectively, using the largest volume data. 



data in the different volumes tend to align. The data with 7 = 0.4 show a good alignment, which 
is taken as the signal of the scaling, eq. The scaling again disappears beyond 7 = 0.4. Thus, 
the optimal 7 is around 7 = 0.4. This analysis is based on the method in Ref. Jlsl ] for an estimate 
of the 7*. 

We attempt to determine the 7^ by fitting the data to the finite-size hyperscaling with the linear 



assumption, eq. (2.3). Since this assumption is valid in large L region, we restrict ourselves to use 
the data in the larger x region. From the fit shown by the dashed line in Fig. |3| we obtain 7^ 0.44 
with ;^^/d.o.f. = 4.0. The value of is reasonably consistent with the value estimated above and 
also that obtained from the hyperscaUng fit in Sect. |]l|. We note that our 7* roughly agrees with 
the result and the results with similar analysis [jl^, |l^] using the data in Ref. 

These results are encouraging as preliminary results, while at present we have not evaluated 
errors of the 7^. There are possibly several systematic errors in this analysis, such as the one due to 
the assumption of the asymptotic form of f{x), corrections of the finite-size hyperscaling from large 



Mf and small L effects [|18[]. If the effect of these corrections is large, we cannot obtain the correct 
even if the data shows the scaling behavior. In order to get reliable estimate of the systematic 
errors, we need to expand our simulation towards lighter mass and larger volume. Detailed analysis 
using those data are ongoing. 





1 1 1 






• L=12 






■ L=18 






• L=24 








♦ 












1^0.0 J 





1 1 1 1 






• L=I2 






■ L=18 






♦ L=24 








♦ 


■ 

1 1 1 1 1 11 . 


-ltO.2 J 



L=12 
L=18 
L=24 



1^0.4 



L=12 
L=I8 
L=24 



-^0.6 



Figure 2: The scaling behavior of Lm^t in the 12 flavors as a function ofx = Lni^^^ where the values of 7 are 
taken as 7 =0, 0.2, 0.4, 0.6 from left to right. The different symbols denote the data at the different volumes. 
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Figure 3: Lm„ in the 12 flavors and the fit result of the finite-size hyperscaling denoted by the dashed line. 
The different symbols denote the data at the different volumes. The filled and open symbols represent the 
data included in the fit and omitted from the fit, respectively. 
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Figure 4: Lnij^ in the 16 flavors at j3 = 3.15(left) and j3 = 3.5(right). The different symbols denote the data 
at the different volumes. The dashed line denotes the fit result of the finite-size hyperscaling. The filled and 
open symbols represent the data included in the fit and omitted from the fit, respectively. 



4. Results of 



/ 



16 



In the Nf = 16 SU(3) gauge theory we perform the simulations at j3 = = 3.15 and 3.5. 
We choose three spatial volumes L = 8, 12, 16 at j8 = 3.15, and add one more volume L = 24 at 
j3 = 3.5. The range of the fermion mass is 0.03 <mf< 0.2. The typical length of the trajectory is 
roughly 1000. 

We analyzed the results of the rrijc and fjc in the Nf = 16 case using the same procedure as in 
the Nf = 12 case discussed in the previous section. The ChPT fit failed as it did for Nf = 12, which 
is expected as the Nf = 16 theory is much deeper in conformal window if Nf = 12 was. 

Using the finite-size hyperscaling analysis, we attempt to determine the value of the at 
two jS's. Figure ffl shows the fit result of the Lnijc with the finite-size hyperscaling assuming the 



asymptotic form of f{x), eq. (2.3). The both panels show good scaling behavior of the Lniji, and 
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the sizes of the ;^^/d.o.f. are reasonable. From the fits we obtain 7* 0.42 and 0.35 at j8 = 3.15 
and 3.5, respectively. We also analyze the data of f^c and obtain 7* which is roughly consistent with 
those with m^j. Again we have not yet estimated errors of the obtained discussed in the above 
section. 

These observations lead to a converged picture of the conformal theory with decreasing 7* 
towards weaker bare coupling. One question arises as to the size of the 7* being much bigger than 
the result with 2-loop perturbation theory y^^^ 0.025. One possible scenario is that the 7* further 
decreases towards much weaker coupling (continuum limit) and eventually gets compatible to the 
perturbation theory. Another possible scenario is our bare gauge coupling is too large to investigate 
the property of this theory in the continuum limit. To investigate further, more detailed study of the 
j8 dependence of the 7^ is necessary. 

5. Summary and outlook 

We have studied the SU(3) gauge theories with the fundamental 12 and 16 fermions using a 
HISQ type staggered fermion action, and presented preliminary results in this report. For the 12- 
flavor case, we attempt to determine the phase of this theory through the analysis of pion mass. Our 
present data favors the conformal hyperscaling over the ChPT The mass anomalous dimension, 7* 
at the infrared fixed point was estimated though the (finite-size) hyperscaling analysis. So far the 
size of 7^ is not as big as 7,, ~ 1 for the theory to be close to the realistic technicolor theory. We have 
not yet estimated the errors of the 7*, so that we need to estimate it to reach a definite conclusion. 
More detailed analyses with more data at larger volume and lighter mass are underway. For the 16- 
flavor case, our data shows the conformal signal, while the obtained value of the 7 is much higher 
than the perturbative result. We plan to investigate this by studying the j3 dependence of this value. 



Acknowledgments 

Numerical calculations for the present work have been carried out on the cluster system "9" 
at KMI, Nagoya University with the public code provided by the MILC Collaboration. We thank 
Katsuya Hasebe for the useful discussion. This work is supported in part by Grants-in-Aid for Sci- 
entific Research from the Ministry of Education, Culture, Sports, Science and Technology (Nos. 
22224003, 23540300, 21540289) and Grants-in-Aid of the Japanese Ministry for Scientific Re- 
search on Innovative Areas (No. 23105708). 



References 



[1] LatKMI Collaboration, K. Nagai et. al.. Many flavor QCD as exploration of the walking behavior 
with the approximation IRflxed point, PoS Lattice 2011 (2012) 080. 

[2] Y. Iwasaki, K. Kanaya, S. Sakai, and T. Yoshie, Quark conflnement and number of flavors in strong 
coupling lattice QCD, Phys. Rev. Lett. 69 (1992) 21-24. 

[3] Y. Iwasaki, K. Kanaya, S. Kaya, S. Sakai, and T. Yoshie, Phase structure of lattice QCD for general 
number of flavors, Phys. Rev D69 (2004) 014507, [[hep-lat 7030 915 9[]. 



6 



Many flavor QCD with Nf = 12 and 16 



Yasumichi Aoki 



X.-Y. Jin and R. D. Mawhinney, Lattice QCD with 8 and 12 degenerate quark flavors, PoS LAT2009 



(2009)049, [arXiv: 0910. 32161. 



Z. Fodor et. al.. Twelve massless flavors and three colors below the conformal window, Phys. Lett. 



B703 (201 1) 348-358, [arXiv: 1104 . 3124| 



A. Deuzeman, M. P. Lombardo, and E. Pallante, On the spectrum of QCD-like theories and the 



conformal window, arXiv: 1201 . 1863 



T. Appelquist, G. T. Fleming, and E. T. Neil, Lattice Study of the Conformal Window in QCD-like 



Theories, Phys. Rev Lett. 100 (2008) 171607, [arXiv: 0712 .0609] 



T. Appelquist, G. T. Fleming, and E. T. Neil, Lattice Study of Conformal Behavior in SU(3) 



Yang-Mills Theories, Phys. Rev D79 (2009) 076010, JarXiv : 90 1 . 3766| ]. 

A. Hasenfratz, Infrared flxed point of the 12-fermion SU(3) gauge model based on 2-lattice MCRG 



matching, arXiv :1106.5293 



T. Aoyama, H. Ikeda, E. Itou, M. Kurachi, C.-J. Lin, et. al, Conformal flxed point ofSU(3) gauge 



theory with 12 fundamental fermions, arXiv:1109.5806 



HPQCD Collaboration, E. Follana et. al.. Highly Improved Staggered Quarks on the Lattice, with 



Applications to Charm Physics, Phys. Rev D75 (2007) 054502, [hep-lat/0 6100 92] 



A. Bazavov, T. Bhattacharya, M. Cheng, C. DeTar, H. Ding, et. al.. The chiral and deconflnement 



aspects of the QCD transition, arXiv :1111.1710 



V. A. Miransky, Dynamics in the conformal window in QCD like theories, Phys. Rev. D59 (1999) 



105003, [ |hep-ph/ 98 123501 1. 

L. Del Debbio and R. Zwicky, Hyperscaling relations in mass-deformed conformal gauge theories. 



Phys. Rev D82 (2010) 014502, [arXiv:1005.23711 



L. Del Debbio, B. Lucini, A. Patella, C. Pica, and A. Rago, Mesonic spectroscopy of Minimal Walking 



Technicolor, Phys. Rev D82 (2010) 014509, [arXiv: 1004 . 3197| 



T. DeGrand, Finite-size scaling tests for spectra in SU(3) lattice gauge theory coupled to 12 
fundamental flavor fermions. 



arXiv: 1109 . 1237 



T. Appelquist, G. T. Fleming, M. F. Lin, E. T. Neil, and D. A. Schaich, Lattice Simulations and 



Infrared Conformality, Phys. Rev D84 (201 1) 054501, [ |arXiv : 1 1 6 . 2 1 4 8| ]. 

Y. Aoki, T. Aoyama, M. Kurachi, T. Maskawa, K.-i. Nagai, et. al. Study of the conformal 



hyperscaling relation through the Schwinger-Dyson equation, arXiv: 1201 . 4157 



7 



